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Abstract 

The purpose of this study is to determine the effect of pressure losses and actual heat transfer on the 
performance of a solar Stirling engine. The model presented includes the effects of both internal and 
external irreversibilities of the cycle. The solar Stirling engine is analyzed using a mathematical model 
based on the first law of thermodynamics for processes with finite speed, with particular attention to the 
energy balance at the receiver. Pressure losses, due to fluid friction internal to the engine and mechanical 
friction between the moving parts, are estimated through extensive and rigorous use of the available 
experimental data. The results of this study show that the real cycle efficiency is approximately half the 
ideal cycle efficiency when the engine is operated at the optimum temperature. © 1999 Elsevier Science 
Ltd. All rights reserved. 
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1. Introduction 

Evaluation of the losses due to the irreversibilities in the Stirling cycle is a topic of 
significant interest to those concerned with the analysis and performance of thermal machines. 
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Nomenclature 

A absorbtivity 

a coefficient 

C concentration ratio 

c v constant volume specific heat 

D diameter 

E insolation 

m mass of the gas 

N number of gauzes of the matrix, number of regenerators per cylinder 
77 r rotation speed 

p pressure 

A p pressure loss 

Q heat transfer rate 

R gas constant 

S surface area 

5 stroke 

T temperature 

AT temperature difference 

u convective heat transfer coefficient 
V volume 

W power output 

w piston speed 

X fraction of heat lost due to incomplete regeneration 
Greek symbols 

e compression ratio, emisivity 

y specific heat ratio 

rj efficiency 

rj' efficiency defined in Eq. ( 6 ) 

p density of the gas 

ct constant of Stefan-Boltzmann 

Subscripts 
aver average 

C cylinder, related to the Carnot cycle 

c related to the working fluid temp. 

cav cavity 

cone concentration 

cpr compression 

exp expansion 

ext external 

f friction 
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g gas 

H source 

L sink 

M material 

mir mirror 

o ambient 

open receiver opening 
R regenerator 

Rec receiver 

SE Stirling engine 
v volumetric 

w related to the piston speed or to the working fluid temperature 

I state in cycle 

II related to the Second Law 


It is the focus of this study. Our studies and recent investigations by others [1,2] have found 
that irreversibilities in the thermodynamic cycle have significant importance in predicting the 
performance of Stirling engines. Efforts have been made over the past several years to better 
understand these irreversibilities related losses [3,4], This has resulted in a number of models 
for the analysis and optimization of Stirling machines that include the effect of irreversibilities 
in the cycle. However, analyses that are based on the conventional entropy or exergy 
techniques do not relate the irreversibilities to the physical phenomena that cause them. The 
model presented here directly connects the irreversibilities to the operation of the cycle at finite 
speed. It provides a clearer understanding of the loss mechanisms and relates them more 
quantitatively to the thermodynamic irreversibility terms. This in turn provides increased 
insight to the loss mechanism. 




Fig. 1. Stirling engine cycle with irreversibilities. 








1726 


M. Costea et at. / Energy Conversion & Management 40 (1999) 1723—1731 


The model presented is an extension of previous work [5,6] and includes the effects of both 
internal and external irreversibilities. Attention in the analysis presented is given to the effects 
of (1) heat transfer through a temperature difference at the source and sink, (2) incomplete 
heat regeneration, (3) piston speed and (4) effects of internal mechanical and fluid friction. This 
presentation extends previous work by basing the analytical pressure loss prediction due to 
internal fluid and mechanical friction on actual operating data and by basing the heat transfer 
rate to the engine on the available heat at the solar source. 


2. Analysis of the solar Stirling engine cycle with irreversibilities 

2.1. The Stirling engine cycle with irreversibilities 

The ideal Stirling engine cycle is shown on p-V and T-S diagrams in Fig. 1. The T S 
diagram is modified to include the effects of heat transfer through a temperature difference at 
the source and sink and for incomplete regeneration. 

Additional heat from the external source is shown to be needed in the process, Q x _ 3 , due to 
incomplete regeneration. Similarly, the unregenerated heat is shown being rejected in the 
process, Q v - \ . Fluid friction of the gas as it passes through the regenerator causes most of the 
pressure loss. 

2.2. Mathematical model 


The solar Stirling engine is analyzed using a mathematical model based on the first law of 
thermodynamics for processes with finite speed, with particular attention to the energy balance 
at the receiver, using what we call “the direct method” [7-10]. 

The net power output of the engine is: 

^ = '?se • Qn — — 7 ( ! — /| \) ^ n ’ irrev ' (1) 


with A T h = T h -T w , and \T L =T C -T L . 

In Eq. (1) the engine efficiency is separated in two terms, the second-law efficiency, Tjn.irrev 
[11] and the Carnot cycle efficiency: 


T l + AT l 
Th-AT h 


( 2 ) 


The heat flux transferred to the working gas in the receiver of the engine has the form: 
Qh = [mR(Tu - AT H ) ■ In e v + X ■ mc v (T H — AT H - Ty — A 71)] • «, 


P\Vi 

t l + at l 


(Th — A 71) 


In e v + 


X 

y -1 


he, AT 


«r 


( 3 ) 


The effect of each of the engine irreversibilities is evaluated independently [10] so the effect 
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of each is clear in the overall optimization procedure. Thus, the second-law efficiency has the 
form: 


'/il.irrev '7ll,irrev int(x) ' ? 7 II,irrev intMp) 


(4) 


where the effect of incomplete regeneration, internal mechanical and fluid friction and piston 
speed are given by: 


>lu 


irrev int(A - } 


1 + 


X 

(}' - 1) • In e, 


’ 1 1c,at 


(5) 


'7 ii 




- 1 


3)1 ■ J2( a p0/pi 




ATu 


T l + AT l 


In Sy 


with 


( 6 ) 


- ’lc.AT ■ '?IT,irrev, 

(7) 

h = 1 , 

3 • £ v 

(8) 

Pressure losses considered in the model are [7-10]: 


^2,(Ap)i = Ap. dverg + Ap f + Apy, 

(9) 

Pressure losses, due to fluid friction internal to the engine were 
and rigorous use of the available experimental data [12]: 

estimated through extensive 

15/1 

d/>aver g — Ap t hrot R — ^ 

(10) 

where the speed of the gas in the regenerator is considered to be proportional to the average 

piston speed: 


w Dp 

UR ~ 2 N R Di 

(11) 


and the gas density in the regenerator, p R , were calculated using the average values of 
temperature and pressure of the gas in the regenerator. 

We note that the pressure losses in the heater and cooler were neglected. Thus, the term of 
Eq. (10) takes account only of the pressure losses due to the throttling of the gas flowing in the 
regenerator. This expression was found by an analytical analysis of the experimental data from 
[12, Fig. 5.7], 

The mean effective pressure losses resulting from mechanical friction of the engine 
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components were estimated. The estimates are based on experimental data from internal 
combustion engine [13] adapted to actual Stirling machine operation: 


Apt = 


(0.94 + 0.045 • w) ■ 10 5 
3)i 



1 

£ v 


[N/m 2 ] 


( 12 ) 


with w [m/s]. 

Pressure losses, due to the piston speed [7-10], actually to the pressure waves generated by 
the finite speed movement of the pistons, waves traveling with the speed of sound in the gas, 
are: 





a • w cpr 

' f'averexp 

^cpr 


^ ' Wexp 

^exp 


where 


c = V3RT- average molecular speed; 


(13) 


(14) 


a = ^J3y —coefficient which depends on the specific heat ratio; w cpi .=H’ ex p = w; /r lvcl - tpr , 
Paver exp — average pressure during the compression and expansion of the gas respectively. 

By combining Eqs. (9)—(14), the pressure losses becomes function of w and T H : 




15 


1 p\ • (1 + £ v ) VV 2 

2R-(T l + AT l ) 'T 


•TV- 


D 2 

NrD 2 Ri 
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(0.94 + 0.045 • w) ■ 10 5 
3)i 


, 1 V 1 £ v 

1 -+ -W ■ P\ -- 

£v / 2 £v — 1 


(15) 


y 


In £ v •, / — • 


1 


R JT l + AT l 


1 + 


Tn-AT U 
T l + AT l 


Taking account of the relation between the piston speed and the rotation speed (s-stroke): 


77 r = 


w 

2s 


(16) 


Eq. (15) becomes a function of n r and T H - 

When Eqs. (3)-(8) and (15) are substituted into Eq. (1), the net power output can be shown 
to be a function only of n r and T H . 

The second assumption of the model is that the heat flux received by the engine is equal to 
the heat flux supplied by the solar receiver: 

Q avail Q H 

where Q H is given by Eq. (3) as a function of the same variables n r and T H . 

The term g aV aii was determinated for a dish mirror with cavity type receiver: 


( 17 ) 
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avail */Rec ' ^cav ’ E ' C ■ )/ c 


(18) 


where the receiver efficiency, r| Rec , depends on T H and optical characteristics of the 
concentration system [5,6,14], 

Finally, Eq. (2) imposes the rotation speed as a constraint in correlating the heat flux to the 
cycle gas with the heat supplied by the receiver. Upon substitution into Eq. (18): 


PiV _i 

Ty + ATy 


■ (T h + AT h ) ■ 


In s v 
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^cav (7 h - To)S C av 
A-cav ‘ J /conc ' E ■ C A cdv ■ )( conc ■ E ■ C ■ Sopen 
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^cxt( Uxt T) fS'cxt 
cav ‘ J /conc ' E ■ C ■ S^pen 


(19) 


where the results have been determined in previous work [5,6,14], 

The analysis to accomplish the optimization of the solar Stirling engine results in a system of 
nonlinear equations that are solved using the method of Lagrangian undetermined multipliers. 


3. Results 

The preceeding analysis of the Stirling engine provides an expression for the optimum speed 
of engine rotation and for the power output, including the effect of irreversibilities, as a 
function of the source temperature. Calculations have been made for illustration using helium 
as the working fluid and with the following fixed parameters: D mirr = 6.5 m, d Re c = 0-4 m, 



Fig. 2. Real power output and the corresponding rotation speed of the engine versus temperature (helium, 
AT = 25 K, ^Rec~ 0.4 m, E = 750 W/m 2 ). 
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r|conc = 0.665, U\ = l W/nrK, u 2 = 6 W/m 2 K, e m = 0.8 and A cav = 0.94, E = 750 W/m 2 , 
T 0 = 293 K, a = 5.67 W/nrK 4 , e v = 5, X = 0.2, j = 0.1m. 

The results of the calculations are shown plotted in Fig. 2. The maximum value of the power 
output IF[fF] of the Stirling cycle engine indicates the optimum temperature of the source (in 
this example /„ [rot/min] is seen to be at about 1000 K). 

The efficiency of the Stirling engine operating with the same fixed parameters, both ideal and 
including the irreversibilities considered in this analysis, is shown in Fig. 3. This figure indicates 
that the ratio of actual to ideal efficiency increases as the temperature increases and that the 
ratio is approximately one half at the optimum temperature for maximum power output. 

The effect of irreversibilities on each of the cited efficiencies as a function of source 
temperature is shown in Fig. 4. The fixed parameters are the same as in the previous 
illustrations. 

It can be seen that for low temperature of the receiver the effect of pressure losses is 
important. The efficiency associated to the incomplete regeneration, r) n irrev int (x), does not 
influence substantially the reversible cycle efficiency in comparison with the pressure losses in 
the range of the optimum temperature we have obtained (Th = 950h- 1100 K, see Fig. 2). This 
conclusion results from the calculation we have made for X = 0.2; 0.4; 0.6. 

The performance predicted (30% efficiency) in the analysis presented is in close agreement 
with actual performance data obtained experimentally from a Stirling engine powered by a 
solar dish [15]. 


4. Conclusion 

An analysis of the solar Stirling engine with irreversibilities is presented that makes possible 
the prediction of the performance of actual solar Stirling engines operating under the same 
conditions. 



Fig. 3. Efficiency of the engine operating upon an ideal or an irreversible cycle versus temperature (helium, 
AT = 25 K, ^Rec “0.4 m ,E= 750 W/m 2 ). " 
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Fig. 4. Efficiencies associated with each class of irreversibility. 
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